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ABSTRACT 

We investigate conditions for and consequences of spallation in radio-quiet Seyfert galaxies. 
The work is motivated by the recent discovery of significant line emission at 5.44 keV in Suzaku 
data from NGC 4051. The energy of the new line suggests an identification as Cri Ka emission, 
however the line is much stronger than would be expected from material with cosmic abundances, 
leading to a suggestion of enhancement owing to nuclear spallation of Fe by low energy cosmic 
rays from the active nucleus. We find that the highest abundance enhancements are likely to take 
place in gas out of the plane of the accretion disk and that timescales for spallation could be as 
short as a few years. The suggestion of a strong nuclear flux of cosmic rays in a radio-quiet Seyfert 
galaxy is of particular interest in light of the recent suggestion from Pierre Auger Observatory 
data that ultra-high-energy cosmic rays may originate in such sources. 

Subject headings: galaxies: active - galaxies: individual: NGC 4051 - galaxies: Seyfert - X-rays: galaxies 



1. Introduction 

Recent X-ray data from XMM-Newton, Chan- 
dra and Suzaku has led to the discovery of nar- 
row line emission in the 5-6 keV regime in sev- 
eral local Seyfert ga laxies (e.g. Turner et al. 20o3 
Yaaoob et al. 20031 ). One popular interpretation 
of the lines has been as emission from so-called 
'hotspots' on the accretion disk, i.e. enhanced 
emission related to events such as magnetic recon- 
nections on the disk surface. The observed line 
profile is modified by Doppler and gravitational 
effects, depending on the inclination of the sys- 
tem to the observers line-of-sight and the radial 
location of the hotspot. A hotspot that co- rotates 
with the disk should show a periodic pattern of 
variations in line energy and strength with time 
(so long as the disk is not observed face-on). If 
the hotspot originates within 20rg then general 



relativistic effects are predicted to be measurable 
using current X-ray data. Superposed on the peri- 
odic shifts, lines should show down-shifting in the 
peak energy of the line with time as the mate- 
rial spirals inward (jPovciak et al. 20041 ). As disk 



hotspots are not expected t o survive for more than 
a few orbits at small radii (jKaras et al. 200 11 ) . ob- 
servation of persistent lines apparently originating 
close to the event horizon would disfavor the disk 
hotspot model. 

Alternatively, the shifted lines may arise in 
ejec ted blobs of gas (s uch as comprising a wind, 
e.g. [Turner et al. 2004 ). In contrast to the hotspot 
model, lines originating from ejecta are expected 
to show non-periodic evolution with time as the 
gas moves outwards, and by tracing the line en- 
ergy over time one can potentially constrain the 
acceleration or deceleration mechanism. In prin- 
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ciple then, one could distinguish disk hotspot and 
ejecta origins for hne emission in the 5-6 keV 
band using time-resolved spectroscopy. 

However, another model has been suggested, 
where observed lines may be identified as species of 
elements such as Cr and Mn that normally would 
be too weak to measure using current data, but 
whose strength has been enhanced owing to abun- 
dance changes in the emitting ga s from spalla- 
tion of Fe (e.g. [Turner et al. The interac- 
tion of protons having kinetic energy > 30MeV 
with matter can result in spallation of its heavy 
nuclei, creating enhanced abundances of elements 
lower in mass than the target nucleus. In the 
low energy regime the protons required for ef- 
fective spallation need be only mildly relativis- 
tic, easily achievable in a number of astrophysical 
situations and therefore we might expect to ob- 
serve spallation under a range of conditions. In- 
deed, spallation is known to significantly affe ct the 
abundance ratios in our own Galaxy ( Reeves 1974 ; 
Lund 19891) where there is approximate energy 



equipartition between cosmic-ray protons and the 
magnetic and radiation fields. The most notice- 
able effects of spallation are that abundant nu- 
clei such as C, N, O, Fe are broken down, increas- 
ing the fraction of lighter elements such that the 
emission or absorption profile of the gas is signif- 
icantly different than expected for cosmic abun- 
dance material. However, the cross-sections are 
relatively low (a < 1000 mb, Letaw et al. 19831: 
Silberberg et al. 19981; iTripathi et al. 1999r where 



1000 mb is 10~^'*cm^), and significant abundance 
changes in a large mass of gas require either high 
proton flux and/or long timescales. 

The conditions for and the consequences of cos- 
mic ray production in AGN have been discussed in 
the past by several authors (Kazanas & Ellison 19& 
[Axford 1994; Cronin 2005 ; Nagano fc Watson 200" 
While the effects of spallation could be detectable 
in any waveband, the cross-section for spalla- 
tion incr eases with atomic mass approximately 
as A'^^'^ ( Letaw et al. 1983 ); this, the high cos- 
mic abundance of Fe and its observability over 
a wide range of ionization mean that signatures 
of spallation may be most easily detected in X- 
ray spectra of iron group elements. At low cos- 
mic ray energies the primary products from the 
spallation of Fe, and their partial cross-sections 
for prompt production following collision with a 



100 MeV proton are Mn (400 mb), Cr (352 mb), 
V (115 mb) and Ti (70 mb), where the cross- 
sections have been determined from code sup- 
plied bv I Silberberg et aT ri998l) . If we adjust the 
cross-sections to account for the decay of unsta- 
ble isotopes, the effective cross-sections become 
Mn (247 mb), Cr (560 mb), V (163 mb) and Ti 
(136 mb), adopting the decay c hannels summa- 
rized in Table 1 of ISkibo ri997h but neglecting 
decay of ^"^Mn, whose half-life is 3.7 x 10^ years. 
At X-ray energies, neutral species produce Ka 
emission at 5.9, 5.4, 4.9 and 4.5 keV respectively. 
Comparison of the partial cross-sections shows 
that the most notable abundance enhancements 
produced from spallation of Fe would be for Cr 
and Mn. 

Spallation of Fe has been discussed previously 
to explain line emission at 5.6 and 6.1 keV as emis- 
sion from ionized species of Cr and Mn in NGC 
3516, although in that case the li nes ratios did not 
agree well with those pr edicted bvlSkibo (19971) for 
spallation of disk gas (ITurner et al. 20021) . How- 
ever, the simulations bv lSkibo fl997l ) were conser- 
vative regarding the level of abundance enhance- 
ment that co uld occur as a result of spallation in 
active nuclei. ISkibo (19971 ) considered the require- 
ments for abundance enhancement in the total 
mass of material surrounding an accreting black 
hole, requiring a high efficiency of proton creation 
to achieve high factors of enhancement for Cr and 
Mn abundance. 

In recent years our understanding of the nu- 
clear environs of Seyfert galaxies has evolved, with 
ever-increasing evidence for large columns of gas 
covering a wide range of ionization state shroud- 
ing the active nucleus and shaping the observed 
properties of AGN in the X-ray bandpass. In this 
paper we argue that, if the gas in the putative 
disk wind provides the target for proton spalla- 
tion, rather than the body of accretion disk itself, 
then the predictions and consequences of the pro- 
ce ss could be v ery different from the calculations 
of lSkibo (199% . 

Recent Suzaku observations of NGC 4051 re- 
veal a line at 5.44 keV that is steady in flux and 
energy over years ([Turner et al. 20101 ) disfavoring 
disk hotspot and ejecta models for the line ori- 
gin and motivating a reconsideration of the im- 
portance of spallation in AGN. Exploration of the 
conditions for and consequences of spallation in 
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AGN is the topic of this paper. 

2. Summary of Observational Results 
from NGC 4051 

Line emission has recently been found at 5.44 
keV and 5.95 keV in 2005 and 2008 Suzaku obser- 
vations of NGC 4051 (iTurner et al. 201Q) . While 
potentially of interest with regard to an iden- 
tification as Mn emission, the line at 5.95 keV 
line suffers a degree of contamination from the 
slightly broadened component of Fe K g emission 
(|Miller et al. 20091 IXurner et al. 2010f ) and (to a 
lesser degree) the detector calibration source. De- 
tailed analysis also showed the inferred signifi- 
cance of the line at 5.95 keV to be sensitive to 
the form of the continuum model. For these rea- 
sons we concentrate on the line at 5.44 keV and 
the inferred ratio of Fe/Cr in this work. 

Principal components analysis can be used to 
decompose data into a mathematical solution 
showing the simplest constant and variable or- 
thogonal vector set that can explain the observa- 
tions. We applied this technique to the combined 
2005 and 2008 data from Suzaku observations of 
NGC 4051. That analysis showed that a line at 
5.44 keV and the neutral component of Fe Ka 
emission comprise part of a flat 'offset' spectral 
component that dominates the sou rce spectrum 
at low flux levels (jMiller et al. 20091 ). 

The statistical case for a line at 5.44 keV is 
very strong in NGC 4051. The improvement to 
the spectral fit using a Gaussian line component 
to model the emission is Ax^ — 32 correspond- 
ing to a detection at > 99.9% confidence; this 
was obtained using a local continuum parame- 
terization and fitting the 2005 Suzaku data. In- 
stead of using a local continuum parameteriza- 
tion one can derive a complex model for the full- 
band Suzaku and HETG data and again test for 
the presence of an additional line. We found 
that such an approach does not diminish the 
level of confidence at which this line is detected 
([Turner et al. 2010t [Lobban et al. 2010l ). The fit 
yields a line energy (in the rest-frame of the 
host galaxy) E = 5.44 ± 0.03 keV, line fiux n = 
5.03^2 01 X 10^^ photons cm~^s~^ [la uncertain- 
ties are quoted throughout) and equivalent width 
(against the total continuum) 46 ± 16 eV for 2005 
Suzaku data. In that fit the flux and equivalent 



width of the Fe Ka line component (observed at 
an energy E — 6.410 ± 0.015 keV) were n = 
1.59±0.23xl0"^photonscm-2s-i and 195±24eV 
respectively. As the width of the 5.44 keV line was 
not determinable with existing data, it was held at 
cr = 50 eV as determined for Fe Ka (whose fitted 
width was a — 50^33 eV) and justified based on 
the evidence from principal components analysis 
that the Fe Ka emission a nd the 5.44 keV em is- 
sion have a common origin (jMiller et al. 20091 ). 

The Suzaku data were time-sliced and the line 
found to be detected in the six resulting time 
periods, conclusively demonstrating that it can- 
not be a statistical fluctuation in the mean spec- 
trum. Further to this, the feature was deter- 
mined to be present in all three operational CCDs 
([Turner et al. 2010f l. 

Finally, we note that the background count 
rates contributes only 2.5% of the total count rate 
in the 5-7 keV band for the summed XIS spec- 
tral data during the 2005 observation when the 
source is relatively dim (and only 1.3% when the 
source is brighter, during 2008). There is no ev- 
idence for line emission at 5.44 keV in the back- 
ground spectrum and the combined background 
spectra from 2005-2008 yielded an upper limit (at 
90% confidence) on the flux of such a line to be 
n < 4.14 X 10"* photons cm"^s~\ i.e. < 1% of 
the detected line flux. The lack of a feature at 
comparable flux or equivalent width in the back- 
ground data also rules out an origin of the 5.44 keV 
line as a detector feature. 

The flux of the line at 5.44 keV was found to be 
consistent with a constant value across the time- 
sliced Suzaku data and in archived XMM data 
( Turner et al. 2010[ ): thus the line shows up most 
prominently when the continuum level is low (Fig- 
ure [T]) , supporting the result from principal com- 
ponents analysis. The line energy was also found 
to be consistent with a constant value and, as 
noted previously, these observations thus disfavor 
hotspot and ejecta models for the line. 

In the rest-frame of the host galaxy the line may 
be identifled as Ka emission from Cr i, in which 
case the observed line strength exceeds what might 
be expected from simple continuum illumination 
of material with cosmic abundance ratios. Strong 
Cr emission that appears linked to the Fe emission 
is n aturally explained by the spallation process. 
The lAnders fc Grevesse (19891 ) abundance ratio is 
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100:1 between Fe:Cr; taking the fluorescence yields 
into account (Ype = 0.347, Ycr = 0.282) we would 
expect line equivalent widths relative to the il- 
luminating continuum to be observed in the ap- 
proximate ratio 123:1 for Fe:Cr from photoionized 
gas; folding in the relative K-shell photoionization 
cross sections modifies the estimate to an approxi- 
mate ratio 113:1 for Fe:Cr. In the absence of spal- 
lation, Cr Kof line emission should be immeasur- 
ably weak for most AGN observed with current 
X-ray instruments. Contrary to expectations, the 
observed equivalent width ratio in the low-state 
data for NGC 4051 is Fe/Cr= 4.24 ± 1.57, i.e. 
showing a significant deviation from the expected 
ratio with a factor of ~ 30 enhancement of Cr Ka 
emission relative to the line ratios expected for 
cosmic-abundance material. The degree of spal- 
lation indicated by the ratio Fe/Cr is much more 
pronounced than observed in the Milky Way, but 
more extreme spallation might be expected as very 
different conditions of cosmic ray fiux and mate- 
rial exposure are thought likely to exist close to an 
active nucleus. 

Fitting a simple Gaussian model to the Fc Ka 
emission in NGC 4051 gave constraints on the 
width of the line that suggested an origin for 
the emission lines at r ~ 0.065 light days, or 
2 X lO^"' cm, with the 90% confidence range be- 
ing 7 X 10" - 3 X 10^5 cm. The strength of the 
Fe line can be used to constr ain the column den- 
sity of the emitting region (|Yaaoob et al. 20091 ) 
yielding Nh ^ lO^^cm"^, and in the toroidal re- 
processor model suggests a global covering fac- 
tor Cg ~ 0.9. Detailed spectroscopy of the 2008 
Chandra HETG and 2005-2008 Suzaku data also 
show a large amount of absorbing gas along the 
line-of-sight, covering a range of column densities, 
ionization states and outflow velocities, support- 
ing the high global covering derived from the line 
strength. Combining the observational constraints 
thus yielded a mass estimate M > 4 x 1O~'*M0 for 
the line-emitting gas with 90% confldence ra nge 
5 X 10-5 < M < 10-1 M0 (iTurner et al. 2010h . 

While the tr ue total pro file of Fe Ka may be 
more complex (j Miller et al. 2009) . the limit esti- 
mated from the simple Gaussian model fit to that 
line provides a useful mass estimate for examina- 
tion of the general feasibility of processes in this 
case study. 



3. Discussion 

With the evidence for a strong anomalous line 
in NGC 4051 and the possible identification as Cr 
emission we explore the possible production of cos- 
mic rays in radio-quiet AGN, the dependence of 
spallation on parameters of the system, and the 
consequences of the process with respect to cur- 
rent observational constraints. 

3.1. Production of cosmic rays in AGN 

The high luminosity observed in the X-ray 
bandpass for AGN is thought to be produced from 
ultraviolet photons originating in the inner accre- 
tion disk that are up-scattered to X-ray energies 
(e.g. iHaardt fc Maraschi 1991 ). The requirement 
for an up-scattering mechanism has led to the 
suggestion that a large flux of hot electrons ex- 
ist in the vicinity of the inner accretion disk and 
consequently an accompanying flux of protons or 
heavier nuclei may be assumed to be present. For 
a steep power-law cosmic ray energy spectrum 
most spallation is caused by low energy protons, 
E < 200 MeV, which are only mildly relativistic 
{P ^ 0-6). Shock acceleration could produce such 
a populatio n of protons, or m echanisms such as 
that of Bl andford fc Pavne (1982 ) , where plasma 
is accelerated by a magnetic fleld tied to the ro- 
tating accretion disk, which might then generate 
protons with those energies. 

3.2. Calculation of abundance enhance- 
ments in the "thick target" limit 



Crosas fc Weisheit (19m and ISkibo (19971) 
have previously considered th e possibility of spal- 
lation in AGN. ISkibo (199"7h considered specif- 
ically the spallation of Fe and its detectability 
in X-ray spectra, and in thi s section we follow 
closely his analysis. ISkibol assumed the case 



where gas presents a thick target to protons: 
i.e. all the cosmic ray energy is absorbed or ra- 
diated and there is no diffusion of cosmic rays 
out of the region. At low proton energies their 
primary energy loss is through Coulomb colli- 
sions, with an energy decay length, expressed in 
units of column density, of appr oximately A n — 
1026(£;/GeV)i-4'55/3 atoms cm-2 (jSkibo 1997h for 
protons of energy E and velocity relative to light 
(3 . Thus zones of gas that are opaque to X-ray 
photons, with Nh > 10^** cm~^, also efficiently ab- 
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sorb cosmic ray protons with E < 100 MeV, and 
for steep incident power-law cosmic ray spectra 
the "thick target" approximation is vahd for col- 
umn densities N > 10^^ cm~^ (the validity of the 
thick target approximation is discussed further in 
Section [3?3)) . Within the target, proton directions 
are efficiently made isotropic by Coulomb scatter- 
ing and magnetic field deflections. The propaga- 
tion of cosmic rays through such a target can be 
treated as a diffusion problem. Those cosmic rays 
then impact heavy nuclei, and the rate per target 
nucleus of spallation of species i into species j is 



Rij = in 



a,j{E)J{E)dE, 



(1) 



where Gij {E) is the spallation partial cross-section 
for that reaction and J{E) is the intensity in cos- 
mic rays of kinetic energy E inside the diffusion 
region. In a thick target, J{E) depends on the cos- 
mic ray injection rate and inversely on the mass M 
of material t hrough which the cosmic rays diffuse 
(|Skibo 19971 ). The expectation number of spalla- 
tion events per target nucleus, (riij), can then be 
written as 



47rM 



■/(-E'min, E^s,^, r, (Ty, Ac, A 



inelastic 



) 

(2) 

where r is the length of time the material is ex- 
posed to cosmic rays, Lcr is the rate of AGN cos- 
mic ray kinetic energy output and is the solid 
angle covered by the target gas. / is a function 
integrated over cosmic ray energy whose kernel is 
a function of the input spectrum, here parame- 
terized by a powerlaw in kinetic energy of index T 
with lower and upper cutoffs -Emin, -E-max, the spal- 
lation cross-section and the path lengths due to 



Coulomb and inelastic collisions, Ajnoiastic- iSkibo 
equated r to the expected time for gas of mass 
M to be accreted into the black hole if the mean 
accretion rate is M: r ~ M/M. He also wrote 
Levi = virisco) Mc^, where r]{rjsco) is the effi- 
ciency of conversion of gravitational energy into 
cosmic ray kinetic energy for accretion of matter 
into the innermost stable orbit, from which equa- 
tion[2] may be written as 



An 



-f{Emin, £^max, T, (Jij, Ac, Ai 



inelastic } 



(3) 



with the mass terms conveniently disappearing. 
After taking into account the network of reac- 
tions, [Skiboj found significant abundance enhance- 
ments were achieved for 77(7-1300) ^ 0.05, assuming 
n = An, and that enhancements of Cr compara- 
ble to those inferre d in sectionp were found for 
vi'i^isco) — 0.1. The lSkibo (19971 ) assumptions re- 
sult in a maximum predicted enhancement by a 
factor 15 for Cr Ka emission over that expected 
from unspallated solar abundance material, while 
depleting Fe by a factor 2. This enhancement cor- 
responds to a predicted ratio for the line equiv- 
alent widths Fe/Cr=3, and the observed ratio is 
Fe/Cr=4.24± 1.57, consistent with this value. 

To achieve consistency with the observed line 
ratio, the necessary value for rj is uncomfortably 
large. However, the value of t] is strongly de- 
pendent on the assumptions made and hence the 
timescale needed for significant spallation; thus we 
proceed to explore whether different assumptions 
may yield consistency with observed data with- 
out the need for very high efficiency or very long 
timescales for material bombardment. 

3.3. Spallation in the "thin target" limit 

A key assumption in the above is that all the 
cosmic ray energy is absorbed or radiated by a 
"thick target" . However, the column densities re- 
quired for this assumption to be valid are N > 
10^^ cm~^ at low cosmic ray energies, so we should 
consider the effect of loss of cosmic rays from the 
bombarded region in targets of lower column den- 
sity. In the thick target limit the rate of spalla- 
tion per target nucleus increases with decreasing 
target column density (or mass), and hence for suf- 
ficiently low target mass the predicted spallation 
rate becomes unphysically large, owing to the ne- 
glect of the diffusion of cosmic rays out of the tar- 
get. To solve for the spallation rate while account- 
ing for diffusion losses would require knowledge of 
the target geometry and magnetic field strength 
and structure. However, we may estimate the col- 
umn density at which the thick-target approxima- 
tion breaks down by comparing the calculated col- 
lision rate per nucleus with that expected in the 
thin-target limit, where a nucleus is exposed to the 
incident flux of cosmic rays. Higher energy cosmic 
rays have lower Coulomb losses and diffuse more 
easily from the target, so we should consider the 
cross-section for collisions that convert Fe into Cr, 
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as we are interested in the validity of the approx- 
imation for the purposes of calculating Cr abun- 
dance enhancement. This approach to estimating 
the limits of validity of the thick target approxi- 
mation is rather crude, and it neglects the possibil- 
ity that a target with a very large internal energy 
in magnetic field could hold an enhanced cosmic 
ray density compared with its surroundings, but 
it does give an indication of the typical column 
density at which the thick target approximation is 
likely to break down. 

Calculation of the partial cross-section for 
the spallation of Fe into Cr is described by 
Silberberg et al. (1998h and we use their code 
for that calculation. However, many of the 
products of an Fe-p c ollision are unstable nu- 
clei, so here we follow Skibo fl997f ) and modify 
the cross-sections to allow for decay into other 
species (see also Section[T|). A significant con- 
tribution to the net Cr production arises from 
decay of 50,52,54]y[j-^ j-^g again neglect decay of 
^^Mn owing to its long half- life). Including this 
cross-section in the integration kernel of equa- 
tion[Tl we find a spallation rate per Fe nucleus of 
i?Fo^Cr =i 2 X 10~S(1024cm-V^H)s"\ assuming 
a cosmic ray spectral index F = 2.4, a low energy 
cut-off lOMeV and a cosmic ray luminosity of 
lO^'^ergs"^ irradiating a target shell at distance 
2 X 10^"* cm. We may compare this with the rate 
expected in the thin target limit, for which we find 
i?Fc->Cr — 3.8 X 10^^ s~^, so we expect the tran- 
sition from the thick to thin target regimes to be 
important at column densities A^h — lO^^'cm"^, 
close to the value we infer for the line-emitting 
gas in NGC4051. For lower column densities it is 
essential to use the more conservative thin target 
limit. 

In terms of the efficiency of spallation, the thick 
target limit ensures a maximum number of spal- 
lations occur per cosmic ray, and hence the max- 
imum amount of enhanced material is produced. 
The total amount of enhanced material would be 
less in the thin target case, but the number of spal- 
lations per target nucleus reaches a maximum, so 
the abundance change in the gas being bombarded 
could be more marked in a target of lower column 
density. We return to consideration of spallation 
efficiencies in section[3?5] 



3.4. Efficiencies, masses and timescales 

A crucial parameter in the calculation of abun- 
dance change is the timescale over which nuclei 
are exposed to cosmic rays. From the above we 
can see that a high value of ) may be ob- 
tained either by a high cosmic ray flux or by 
a long timescale, greater than M/M. Assum- 
i ng a bolometric luminosit y L^oi = 10**'^ergs~^ 
(IVasudevan fc Fabian 20091 corrected to the Tully- 

m 



Fisher distance 15.2 Mpc, iRussell 20041 ) and ac- 



cretion bolometric radiative efficiency rf''^^ = 
0.05 we deduce a mass accretion rate M ~ 
0.0035 M0year~^. The mass of gas bombarded 
by the cosmic rays is model-dependent, as we 
discuss below, but for a generic mass of M ~ 
0.01 Mq the accretion timescale M/M is only 
about 3 years. This timescale is much less than 
the current growth timescale of the black hole, 
tbh — Mbh/M ~ 5 X 10^ years for black hole 
mass Mbh = 1-7 x IQ^Mq (|Dennev et al. 20091 ) 
and leaves open the possibility that gas could 
achieve significant enhancements in abundance of 
Cr with a low cosmic ray luminosity, if it can 
survive near the black hole without being either 
accreted or blown out. 

3.5. Location of the enhanced gas 

Arguably the first location to consider for the 
target gas is the accretion disk. However, sim- 
ple models of spallation within the disk both re- 
quire the spallation timescale to be sufficiently 
long to generate the inferred abundance changes, 
T ~ Mdisk/-^, and predominantly generate en- 
hanced abundances only in the inner regions of 
the disk (see Appendix). Spallation may well oc- 
cur within the accretion disk, but it seems un- 
likely that the efficiency would be high enough to 
achieve the abundance changes that we observe, 
and we would require a mechanism that transports 
enhanced material from near the ISCO to larger 
radii for consistency with the line widths that we 
observe. While transportation of the spallated 
gas is possible via outflow mechanisms such as a 
disk wind, the spallated gas detected in NGC 4051 
yields a limit on the outflow bulk velocity v < 1200 
km s~^. It is difficult to determine from X-ray 
spectroscopy whether the component of outflow- 
ing gas traced by He-like and H-like Fe absorption 
lines has been modified by spallation, as Cr would 
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be highly ionized in such a zone and would provide 
little X-ray opacity. 

A more attractive possibility is that the tar- 
get gas is not part of the accretion disk itself, but 
rather is material lying out of the plane of the disk, 
perhaps even in or near the broad line region. In 
this case a relatively small mass of gas could in- 
tercept a large fraction of the cosmic rays emitted 
from the AGN, leading to much higher spallation 
rates than would be inferred for spallation within 
the accretion disk. Detailed timing analysis of the 
NGC 4051 data shows frequency-dependent time 
lags between the hard and soft X-ray photons of 
up to 970s. These may be explained by the effect 
of reverberation in the hard X-ray band as con- 
tinuum photons are reflected from a thick shell of 
circumnuclear material extending to ~ 1.5 x 10^"* 
cm from the black hole, having global covering fac- 
tor > 0.44 ([Miller et al. 200% . This sheU could 
comprise at least some fraction of the target gas. 

To estimate the expected spallation timescale 
in the thick target approximation, we should as- 
sume a high column density, N > 10^^ cm~^, and 
equivalently a relatively high mass of target gas, 
■^gas ^ 4 X lO^^'Mo, for the nominal radius of 
2 X 10^** cm and high filling factor (section[2l re- 
call that in the thick target approximation, the 
abundance enhancement for a given cosmic ray 
exposur e time decre ases with M). Comparing 
with the lSkibo (19971 ) calculation, the exposure to 
cosmic rays would need to last t > Afgas/M > 
4 X 10-4/0.0035 ~ 0.11 years for t] = 0.1° (equiva- 
lent to a cosmic ray luminosity 2 x lO'*'^ ergs"^). 
If we adopt a less extreme value for the efficiency 
of cosmic ray production, we would still only re- 
quire the gas to be exposed to cosmic rays for a 
time T > 1.1/77 years. The reason for the dramatic 
change in timescale, and the less restrictive re- 
quirement for the efficiency rj, compared w ith the 
accretion disk case or with the Skibo (19971 ) analy- 
sis, is the higher rate of cosmic ray bombardment 
per Fe atom in the target gas. A small fraction 
of the gas in the accreting system receives a dis- 
proportionately high fraction of the cosmic rays. 
This of course is quite natural in the standard ac- 
cretion disk picture, where most of the mass is in 
the disk but where additional traces of material 
are exposed to the central source over the remain- 
der of the 47rsr. 

For lower column densities and target masses. 



the timescale needed for significant spallation is 
expected to asymptotically approach the inverse 
spallation rate r — > 1/i? for exposure to the in- 
cident cosmic ray flux (section !?^ . which for Fe 
with the assumed cosmic ray spectrum, luminosity 
lO''^^ erg s"^ and radius 2 x 10^^ cm is t > 0.9 years. 
This estimate for low column densities yields a 
similar value to the thick-target calculation, al- 
beit not one that takes into account the network 



of reactions as in the calculation of lSkibo (19971 ). 

If we knew the detailed geometry and magnetic 
field structure of the target we should solve the 
diffusion-loss equation for some assumed system 
mass and diffusion loss coefficient, but in either the 
thin- or thick-target case the required timescale is 
only of order years, and over the lifetime of the 
black hole growth, a substantial total amount of 
material could be processed through the cosmic 
ray region: at any one instant of time we only 
need to be observing a small mass of spallation- 
enhanced gas which could have been substantially 
enhanced in its Cr abundance. This scenario is 
an effective way of explaining the apparent abun- 
dance enhancements that we see without requiring 
extreme conditions of cosmic ray production. 

3.6. Gamma-ray and radio emission 

Interactions of cosmic rays with the accret- 
ing gas ar e expected to prod uce secondary 7-ray 
emission (iDermer 19861) and ISkibo (19971 equa- 
tions 11-12) estimated the 7-ray flux above 100 
MeV for Seyfert galaxies if spallation is occur- 
ring. We have recalculated the expe cted 7-ray flu x 
specifically for NGC 4051 following [Skibo (19971 ). 
The estimate considers only production of 7-rays 
following neutral pion creation, as calculated for 
the Milky Way bv lDermer (19861 ). whose code we 
use for the c alculation, m odified to incorporate 
equation 4 of ISkibo (1997 ). The calculation as- 
sumes the "thick target" case where all protons 
are absorbed, and hence this provides a maxi- 
mum 7-ray luminosity for a given output AGN 
proton luminosity. Assuming that Lcn = ifeoi, 
the predicted 7-ray flux is i^(> 100 MeV) ~ 3.3 x 
10-®(ri/47r)photons s""'^ cm"^ for proton spectral 
index a = 2.4 and proton low-energy cut-off 
10 MeV. About 0.8 percent of the input cosmic ray 
energy is radiated as 7-rays. For this proton in- 
put spectrum, the 7-ray spectrum is harder, with 
an approximately power-law form with photon in- 
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dex 2.1 in the range 100MeV< < 100 GeV. 

discuss the detection Umit of the 
Fermi Large Area Telescope and for {I ~ 0.9 the 
predicted flux of NGC 4051 falls just below the 10a 
detection limit from the 3-month sky survey data 
~ 3.5 X 10~®photons s~^ cm~^ for a source of this 
spectrum and Galactic coordinates |Abdo 2009^ . 
A "thinner" target in which protons escape the 
AGN, lower global covering factors, lower cosmic 
ray luminosities achieved through longer spalla- 
tion timescales or the possibility of intermittent 
cosmic ray output all weaken the possible con- 
straints available from measurements of 7-ray flux. 

Radio synchrotron emission might also be ex- 
pected, partly as a result of secondary electron 
creation but more significantly through primary 
cosmic ray electrons generated as part of the pro- 
cess that accelerates the protons. Any estimate 
of the synchrotron emission is extremely uncer- 
tain, as we don't know the efficiency of electron 
acceleration compared with proton acceleration, 
we don't know the magnetic field strength and 
the electrons themselves are expected to suffer 
significant ionization losses propagating through 
the material as well as adiabatic and radiative 
losses through bremsstrahlung, synchrotron and 
inverse Compton scattering. Furthermore, any 
compact synchrotron-emitting region would be 
self-absorbed, significantly reducing the detect- 
ability of any synchrotron emission. We do not 
here solve the diffusion-loss equation for electrons, 
but we can already place a stringent limit from 
simple estimates of the radiative losses. The ra- 
tio of synchrotron to inverse Compton losses is 
given approximately by the ratio of energy den- 
sities in magnetic fields, Umag, and radiation, u^^ 
so if these were the dominant energy loss mecha- 
nisms, we would expect an integrated synchrotron 
luminosity 



-'^mag 



neglecting synchrotron self-absorption for the mo- 
ment. If other mechanisms lead to comparable en- 
ergy losses, the synchrotron luminosity would be 
lower. NGC 4051 has an steep spectrum radio core 
with flux density 3.2 mJy at 5 GHz when measure d 
at resolution l.larcsec (|Ho fc Ulvestad 20011 ) . 
corresponding to a spatial resolution 80 pc. The 
core has been claimed to be marginally detected at 



EVN resolution ijCiroletti fc Panessa 20"09l ). The 
core radio luminosity of NGC 4051 integrated up 
to 100 GHz is approximately 3 x lO'^^ergs"^, so 
we can place an approximate upper limit on the 
magnetic field strength if we require Lgync to 
be less than this value. For Leiectron — ^boi — 
lO^^ergs"^ and an energy density in radiation 



at r — 2 X 10 cm of 



700 erg cm * we 



find 



^mag 



< 0.002 erg cm , corresponding to a 



magnetic fiux density B < 0.2 Gauss. In the 
synchrotron-optically-thin regime a discrepant ra- 
dio flux could only be produced if B were higher 
than this value, but then any mJy synchrotron 
source with magnetic field higher than this value 
would be strongly self-absorbed at GHz frequen- 
cies on these angular scales. So even without con- 
sideration of the other energy loss mechanisms, we 
conclude that no significant excess of radio emis- 
sion would be expected from the region. If some 
fraction of relativistic electrons escape to larger 
radii where inverse Compton losses are less se- 
vere and where the synchrotron surface brightness 
limit is not so restrictive, detectable radio emis- 
sion could be produced, but without a detailed 
model for the electron energy losses and electron 
transport to larger radii, we cannot estimate its 
possible contribution to the observed radio emis- 
sion. 

3.7. Energy and momentum transfer 

At low cosmic ray energies < 1 GeV, the cosmic 
ray energy is absorbed by either Coulomb inter- 
actions or inelastic scattering and spallation, and 
ultimately that energy ends up as either heat o r 
kinetic energy (see also Crosas fc Weisheit 19961) . 
For a "thin" target the heating may be rela- 
tively low, but for a thick target we can es- 
timate the energy density in the diffusing cos- 
mic rays by evaluating the cosmic r ay intensity 
J{E) from equation 4 of lSkibo (19971 ). Assuming 
LcR — lO^'^ergs"^ we estimate an energy density 
UcR — 240/A^25 ergcm"^, where N25 is the column 
density in units of lO^^cm"^, assuming a powerlaw 
cosmic ray spectrum with index 2.4 above lOMeV 
and a source-gas distance of 2 x 10^^ cm. The ex- 
pression is invalid for N25 < 1 as the gas would no 
longer satisfy the "thick target" assumption. For 
A^25 = 1 the energy density is equivalent to the 
thermal pressure in gas of temperature 10^ K with 
density 10^^ cm~^, so in sharing their energy with 
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the absorbing gas, the cosmic rays make a signifi- 
cant contribution to the thermal history of a thick 
target. The momentum transfer is similarly sig- 
nificant: in an equivalent of the "Eddington lumi- 
nosity" calculation, for Lcr = lO'^^ ergs~^ and col- 
umn density = 10^^ cm~^ the rate of cosmic ray 
momentum transfer is 0.4 of the Eddington rate. 
These high values of energy and momentum depo- 
sition do not apply in the "thin target" regime, but 
if there were significant column densities of mate- 
rial intercepting the cosmic rays, they would be ex- 
pected to experience si gnificant heating and rad ial 
forces, as suggested bv lSironi fc Socrates 

3.8. On the origin of UHECRs 

The astrophysical origin of cosmic rays is a 
topic of great interest. Cosmic rays up to ener- 
gies of 10^^ eV appear to originate from supernova 
remnants, those in the range 10^^ < E < 10^^ eV 
also appear to have a Galactic origin, however, the 
source of ultra-high energy cosmic rays (UHECR) 
E > 10^^ eV is not yet accounted for. The extreme 
conditions in the nuclei of active galaxies offer an 
appealing possibility for production of UHECRs 
but cosmic rays have not yet been directly estab- 
lished as originating from active nuclei. While the 
cosmic rays required for spallation are of low en- 
srgyj ^ 1 GeV, the mechanisms thought to be in 
play in AGN may also produce UHECRs. 

To date, radio loud AGN have been the focus of 
studies of cosmic-ray production in AGN because 
radio jets have been widely thought to provide a 
suitable site for particle acceleration. However, 
even the jets of radio-loud AGN have fine-tuned 
requirements for successf ul particle acceleration: 
Farrar fc Gruzinov (20091 ) argue that shock accel- 
eration in a relativistic jet is constrained by the 
need to have a sufhciently strong magnetic field 
for acceleration of particles to energies above 10^^ 
eV, while staying within a range that avoids over- 
production of synchrotron radiation and excessive 
photo-pion energy losses in the AGN. 

It had been thought that UHECRs could not 
be produced in the nuclear regions of radio- 
quiet AGN because of the small size of the 
acceleration region and the energy losses (e.g. 
Norman et al. 19"95h . However new work by 
Pe'er et al. (200i^ itias argued that nearby radio- 
quiet AGN could indeed be the source of UHE- 
CRs. Pe'er et al. C2009f ) consider particle accel- 



eration in the parsec-scale weak jets that are 
known to exist in m any local radio-quiet AGN 
including NGC 4051 (|Giroletti fc Panessa 200^ . 
However, recent AUGER results have suggested 
that UHECRs may be dominated by h eavy nu- 



clei (jUnger et al. 20071: iBeUido et al. 200 9). and 



in this ca se the constraints on jet luminosity are 
Pe'er et al. r2009f ) calculate conditions 



reduced. 

for a bolometric luminosity lO^'^ergs"^, applica- 
ble to NGC 4051, and find that the nucleus can 
survive photo-disintegration if the acceleration oc- 
curs on a parsec scale, concluding that radio-quiet 
AGN are viable sources of UHECRs. 

The relationship between the observed photon 
fiux and currently observed cosmic ray flux from 
a particular galaxy is unclear, owing to the signif- 
icant time delay that must exist for the arrival at 
Earth of cosmic rays co mpared with X-ray pho - 



tons (of order 10^ years; Moskalenko et al. 2009f l. 



However one can attempt to conduct a statis- 
tical survey of the coincidence of UHECR oc- 
currence and the positions of AGN. Intriguingly, 
Zaw et al. (200^ have discovered a significant an- 



gular correlation between low luminosity AGN and 
high-energy cosmic rays in data from the Pierre 
Auger Observatory. The correlation is too strong 
to be simply explained by AGN tracing the large- 
scale distribution of matter, indicating that a sig- 
ni ficant fraction of UHECRs are produced in AGN 
(|Farrar et al. 20091). 

Moskalenko et al. r2009l) suggest that as the 



UHECR candidate AGN have no special char- 
acteristics to set them apart from the AGN 
population (i.e. they are not the radio-loud 
sources) then the apparent correlation of low lu- 
minosity AGN with UHECRs m ust be a chance 
occurrence. Zaw et al. (20091) note that the 
candidate AGN may not appear 'special' at 
any given time if cosmic ray acceleration oc- 
curs via a sporadic event such as tidal disrup- 
tion, which may inject > 10^^ ergs int o the 
nuclear environs ( Farrar fc Gruzinov 2009[) and 
yet not be occurring at the time of observa- 
tion of the candidate AGN. Tidal disruption 
events should only occur eve ry 10^ — 10^ years 
( Maeorrian fc Tremaine 1999[ ) in a given AGN 
and the duration of the event itself would likely be 
very brief, of the order of tens of rotations at the 
innermost stable orbit, i.e. hundreds to thousands 
of seconds in the case of NGC 4051. 
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However, current data cannot rule out the pos- 
sibility that modest luminosity, radio quiet AGN 
commonly and persistently produce cosmic rays in 
their nuclear environs. As the typical number of 
UHECR detected is only ~ 1 per AGN, Poisson 
sampling guarantees that not all AGN will have as- 
sociated UHECR in the Pierre Auger Observatory 
data. If UHECR production is linked to produc- 
tion of low energy cosmic rays in AGN, signatures 
of spallation may be better indicators of cosmic 
ray acceleration rather than large-scale radio jets, 
and invoking a sporadic event for cosmic-ray ac- 
celeration appears unnecessary. If spallation is oc- 
curring in NGC 4051, then this provides important 
observational evidence for the existence of a signif- 
icant flux of low energy cosmic-rays in radio-quiet 
active galactic nuclei, and by extension is support- 
ive of the production of UHECRs from the same 
sources. 

4. Conclusions 

Significant line emission at 5.44 keV in NGC 
4051 may be interpreted as evidence for spalla- 
tion at work in the nucleus of a radio-quiet ac- 
tive galaxy. Examining conditions for and con- 
sequences of spallation in circumnuclear material 
we find the highest abundance enhancements are 
likely to take place in target gas away from the 
plane of the accretion disk with timescales for spal- 
lation that could be as short as a few years if the 
cosmic ray output is comparable to the bolometric 
output of the active galaxy. Spallation timescales 
would be proportionally longer for lower cosmic 
ray luminosity, but even so could still be orders 
of magnitude shorter than the accretion timescale 
onto the black hole. Feedback of material from 
the active nucleus to the host galaxy may result 
in enhanced abundance ratios at larger radii than 
probed with these X-ray observations. 

As the acceleration processes required to pro- 
duce the protons necessary for spallation likely 
produces a broad spectrum of energetic particles, 
the results are also interesting in the context of our 
broader understanding of the origin of cosmic rays. 
The evidence for spallation in this radio-quiet ac- 
tive galaxy may support the discovery from the 
Pierre Auger Observatory data that the origin of 
ultra-high-energy cosmic rays may be the nuclei of 
radio-quiet AGN rather than luminous radio-loud 



AGN. 
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A. Appendix material 

A.l. Accretion disk spallation 

We might suppose that material in the accretion disk is exposed to cosmic rays. However, in the simplest 
pictures, spallation wo uld occur pr imarily close to the black hole and would be limited by the accretion 
timescale assumed by ISkibo (199T\) . t ~ M/M. Consider accretion disk material being illuminated by 



cosmic rays from a source a height d above the disk on the rotation axis. An annulus on the disk containing 
mass dM and subtending a solid angle df2 at d, has a cosmic ray intensity 

L dn 

J (X 

An dM 

if the cosmic rays are localized to that radius, as expected for the very high column density in the accretion 
disk. If mass flows through that annulus at a radially-invariant rate M the expectation value for the number 
of spallations per nucleus is 

d{nij) (X 

' 47r M 

and so the total number of spallations per nucleus integrated from infinity to radius r is 

T rl 

dn(r!) - ^ . 

! 2M + r2 

Thus half the expected spallati on events oc cur within a radius ri/2 = VSc?- The total number of spallations 
is half the value calculated by lSkibo (1997 ) owing to the global covering factor of the disk of 0.5, but the 



spallation timescale has the value he assumed. 

We reach a similar conclusion if we suppose instead that cosmic ray sources, perhaps associated with 
shocks in the accretion disk, are embedded within the disk. We gain a factor 2 if the cosmic rays are now all 
intercepted by thick disk material. If the cosmic ray intensity as a function of radius is given by the gradient 
in potential energy, then for Keplerian orbits in a Newtonian potential (j) we expect 

ec^ d(j) 
2 dt 

where e is the efficiency of converting gravitationa l potential energy into cosmic ray energjQ. Note that the 



mass of an element of gas cancels as in lSkibo 19971 The expected number of spallations per nucleus at radius 
r is 

mini -^max7 (^ij i -^inelastic ) (Al) 

which is the relation given by equation[3] for r ~ rigco and = At:. The value of •q{r) should be calculated 
appropriately for some assumed spin of the black hole, but for any metric the spallation occurs predominantly 
close to the ISCO. 



^The relationship between e and 7y{r) for accretion in to radius r is 



2rc2 

assuming Keplerian orbits in Newtonian gravity and 

(1 - 2GMBH/rc^) 



ri{r) 



1 - 



Vl - 3GMBjj/rc2 



for circular orbits in a Schwarzschild metric. Note that the energy released from accretion must be partitioned between all forms 
of energy output from the active nucleus: i.e. it is not possible to have a Schwarzschild black hole producing electromagnetic 
radiation and cosmic rays with ri{risco) = 0.057 in both outputs simultaneously. 
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Thus in either case, we reach the follow conclusions: (i) the spallation timescale is about the same as 
the accretion timescale, r ~ M/M , so to achieve large abundance changes a high efficiency of cosmic ray 
production is needed, rj{risco) ^0.1; (ii) spallation in the accretion disk would be expected to occur primarily 
in the inner regions unless the cosmic ray source is very extended and/or distant from the disk. Yet the Fe Ka 
line, which we suppose originates in the same enhanced gas, is narrow, with a likely location r > 7 x 10^'^ cm 
(section[2]) . However, we don't necessarily observe the gas in the same location as where the spallation 
occurred: thus an accretion disk origin for the spallation would require gas that has been enhanced near the 
black hole to have been transported back out to the large radii where we observe it. 
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Fig. 1. — Suzaku XISO + 2 + 3 data and residuals 
from the mean 2005 spectrum, compared to an 
absorbed powerlaw plus Gaussian line at 6.4 keV, 
showing the residual excess counts at ^ 5.44 and 
5.95 keV 
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